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Abstract 

Data on e + e" — > ir + ir~T(lS,2S,3S) show a large increase in branching fractions near 
T( 10860). A suggestion of Ali et al. [1] is to interpret this as evidence for a tetraquark, 
1^(10890) = (bn)(bn). However, it may also be interpreted in terms of T(10860) —> BB* , 
B*B* and B S B* above the open-6 threshold, followed by de-excitation processes such as 
' BB* — > T(1S, 2S, 3S). In the charm sector, a hypothesis open to experimental test is that 

m X, Y and Z peaks in the mass range 3872 to 3945 MeV may all be due to regular 3 Pi and 

3 P2 cc states (and perhaps 3 Po an d mixed with meson- meson. 

^ ■ PACS numbers: 11.80.Et, , 13.25. Es, 14.40.-n 

S 1 - 1 Introduction 

a,; 

There is much speculation about the existence of tetraquarks. Most papers unfortunately do 
not distinguish between tetraquarks and simple meson-meson. It is prudent, but less exciting, 
to examine conventional interpretations in terms of non-exotic fin, ss, cc or bb states with 
admixtures of meson-meson. 

Nuclei may be viewed as six quark combinations. However, nucleon-nucleon phase shifts 
^ ; may be fitted in terms of meson exchanges. The essential reason nuclear matter does not 
collapse is short-range repulsion due to the Pauli principle. Returning to meson spectroscopy, 
the conventional qq states appear as octets and singlets (or linear combinations). Jaffe has 
suggested that two coloured quarks can form a colour 3 C which is anti-symmetric, or 6 C which 
is symmetric [2]. Then + mesons can be formed as (33) combinations of flavour and colour. 
This neatly evades repulsion between qq and also accounts for the high masses of a (980) and 
/o(980). However, Jaffe's scheme without meson- meson disagrees with observed branching ratios 
^ ' [3]. The ratio g 2 (/o(980) — > K K) / g 2 (ao(980) —> KK) (where g are coupling constants) requires 
that /o(980) has a dominant KK component; also the predicted ratio g 2 (n — > Kir)/g 2 (a — > tttt) 
is much too small. What may defeat Jaffe's proposal is chiral symmetry breaking. The pion and 
kaon are abnormally light, with the result that meson-meson configurations {qq){qq) dominate 
over (qq)(qq) in a, k, a (980) and /o(980). Present measurements of branching ratios are not 
good enough to rule out some small admixture of {qq){qq)', further improvements in the branching 
ratios of a and /o(980) to KK, rji] and tctc above 1 GeV would help greatly. 

Meson exchanges do make good predictions for titi [1] and Kn [5] phase shifts up to masses 
where qq resonances appear. They also predict correctly low energy 1 = 2 and 1 = 3/2 phase 
shifts. There is no evidence for the 27 representation predicted for (6, 6) combinations of flavour 
and colour. In a valuable review, Richard [6] discusses the issues in terms of the flux-tubes within 
four-quark configurations. 

There is a further feature relevant to meson-meson configurations. The train of argument is 
as follows. Many mesons are observed at or very close to thresholds of opening channels. Well 
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known examples are /o(980) and ao(980) at the KK threshold, /2(1565) ar the ujuj threshold 
and X(3872) at the DqD* threshold within ~ 0.3 MeV. The mechanism for this synchronisation 
is rather fundamental [TJ. The conventional form for the denominator of a resonance is 

D(s) = M 2 -s-tY,G 2 tPl (s), (1) 

i 

where s is Lorentz invariant mass squared, Gi = giFi(s), gi are coupling constants of open 
channels and Fi are form factors. However, the correct form for D(s) is M 2 — s — Rell(s) — 
zlmll(s) where Imll(s) = G%(s)pi(s) and 

Ren ( .)=jp r dS 'j:9Mm. (2) 

7T Js thr j S — S 

Here Sthr is the value of s at the opening of the threshold and P denotes the principal value 
integral; [ImII(s) is the pole term from this expression.] The origin of Eq. (2) is that amplitudes 
are analytic functions of s, so that any change in the imaginary part of the amplitude must be 
accompanied by a change in the real part, or vice versa. 

Fig. 1 sketches the behaviour of Imll(s) and Ren(s) at an S-wave threshold, KK in this 
example. The value of Re n(s) is large and peaks exactly at the .S-wave threshold. It acts as 
an attractor and can explain why /o(980) and ao(980) lie very close to the KK threshold. The 
present limitation in calculating Ren(s) accurately is that form factors i^(s) are poorly known. 
Incidentally, Eq. (2) is equivalent to solving Schrodinger type equations, a procedure adopted 
by many authors. These solutions are explicitly analytic. It is also included in the model of van 
Beveren and Rupp, because their amplitudes are constructed algebraically to be analytic [8]. 




0.5 1 

M (GeV) 

Figure 1: Re (s) and 9kkPkk(s) for / (980), normalised to 1 at the peak of g 2 KK PKK- 



The Hamiltonian for a qq state decaying to meson-meson obeys 



H» = ( f l_ ) *; (3) 
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Hn describes short-range qq components and if 22 refers to ingoing and outgoing mesonic chan- 
nels and must include t- and w-channel meson exchanges; V accounts for the coupling between 
them due to s-channel decays. The eigenfunction \1/ becomes a linear combination of qq and 
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meson- meson. The latter is not an 'optional extra'; if a resonance decays to meson- meson, that 
component is a necessary part of the wave function, and indeed plays a vital role. This is the 
basis of a large number of papers by Oset and collaborators, enumerating the attractive meson- 
meson contributions to a large number of mesons, e.g. [9J; their calculations include empirical 
short-range terms which may well simulate qq contributions. 

The form of Eq. (3) is strictly analogous to the formation of covalent bonds in chemistry 
[TO] . According to the variational principle, the eigenstate minimises the eigenvalue. The qq 
component is of short range. Mixing with meson-meson components at longer range lowers 
momentum components in the wave function and hence the energy eigenvalue. 

The J/*&, ^'(IS), T(1S) and T(2S) are very narrow and it is a good approximation to view 
them as pure cc states. Some commentators then argue that other mesons cannot be cc or bb 
because they have unexpected meson-meson components. That is wrong. If a resonance is close 
to an opening threshold with the same quantum numbers, it must contain virtual components 
of that channel. Taking ccnn as an example, there are attractive long range mesonic compo- 
nents in the configurations (cn)(cn) and (cc)(nn); examples are DD* and ojj/^f, which appear 
prominently in observed decays of X,Y,Z in the mass range 3872 to 3945 MeV. The observed 
meson-meson decays are in fact a signature of mesonic components in the wave function. Wein- 
stein and Isgur modelled the /o(980) and ao(980) in terms of meson exchanges [TT] . as did the 
Julich group of Janssen et al. [T2] . 

These general remarks set the scene for further details. Section 2 discusses T( 10860) in 
this light. It might be the first clear tetraquark. However, it may also be understood at least 
qualitatively as a bb state coupled strongly to BB*, B*B* and B* S B* S channels which open in the 
mass range 10559 to 10826 MeV and can de-excite to 7iirT(nS), n < 5. 

Section 3 discusses states popularly known as X, Y, Z. The Particle Data Group [13] labels 
them as x c2 (2P), X(3940) and X(3945); the last of these is somewhat confusing in view of 
the fact that its average mass is now 3915 MeV. Here Z(3930) will be used as a shorthand for 
X c2 (2P). All these states lie close to the DD* threshold. This leads to my suggestion that all 
of these states are n = 2 cc P-states mixed with meson-meson. 

Spin-parity determinations of these states are presently lacking. This appears to be due to the 
fact that only decays have been analysed; that leads to ambiguities. It is essential to analyse the 
full process of production and decay, as in the analysis of Dalitz plots. This leads to orthogonal 
matrix elements between all J p . The ideal process for analysis is B — > K(ojJ/^f); here B and K 
both have spin 0, making the formulae rather simple. Formulae for this process were published 
in 2005 [14] , but experimental groups have not taken advantage of spin information from decays 
of u> and J/^; that is essential because it determines spin alignments. With the full formulae, 
my simulation (which assumes uniform geometric acceptance) suggests that ~ 50 events are 
sufficient to determine J p if only a single amplitude is present or dominant. If there are more 
amplitudes, it gets harder, but it would still be valuable to know the outcome. 

Section 4 reviews the Partial Wave Analysis using tensor notation in its simplest form. The 
formulae can be used in fully relativistic form. However, for the mass range 3872 to 3945 MeV, 
relativistic corrections are small and barely significant with present statistics. Two appendices 
discuss the tensor notation, also the simplest choice of axes and questions of how to do Lorentz 
transformations and rotations of axes, if needed. Section 6 summarises conclusions briefly. 
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2 The T(10860) 

The Belle Collaboration reports branching fractions for T(10860) — > 7t + tt~T(IS), tt + tt~T(2S) 
and 7r + 7r~Y(3S') a factor ~ 50 larger than for T(4S) [15] [16J. Ali, Hambrook, Ahmed and 
Aslam interpret these branching fractions in terms of a tetraquark 1^(10890) = (bn)(bn) [lj. 
More recently, Ali, Hambrock and Mishima [17] predict for this configuration an intensity ratio 

a[T(lS)K + K-]/a[T(lS)K°K } = 1/4. (4) 

For a bb or BB* composition, this ratio should be 1, allowing a clear experimental test. 

There is a natural explanation of the sizable branching fraction for T(10860) — > ttt[T(1S, 2S, 3S) 
in terms of diagrams like that in Fig. 2. The intermediate states are off-shell, allowing n pro- 
duction, even though B* and B* are below the kinematic threshold for pionic decays. The final 
step allows de-excitation of B*B to lower lying 7T7tT channels with larger phase space. Many 
similar diagrams exist with intermediate B, B*, B s and B*. Via mixing with meson-meson 
configurations, the wave function of T( 10860) becomes 

|T(10860) >= a\bb > + Y.( 3 ij\B i B 3 > + £ 7„,|7r7rT(nS) > . (5) 

ij n=l,3 

However, it is not possible to calculate these branching fractions accurately for lack of experi- 
mental information. The decay width of T(4S) to 7T7rT(lS') is very weak, probably because it 
decays almost entirely to BB. Transitions between B and Bit are forbidden by angular momen- 
tum conservation. Channels involving real B* and B* are closed, so only virtual intermediate 
states contribute. However, once intermediate BB*, B*B*, B S B S and B S B* channels are open, 
there are 20 alternative diagrams. There may be complicated interferences amongst them. The 
conclusion is that the observed 7t + tt~T(1S) decays of T(10860) are not dramatically out of line 
with expectation, bearing in mind the increase in phase space. What would be decisive is an 
accurate determination of the branching ratio of T(10860) to K + K~ and K°K°. 

A technical detail concerns the fit to the 7t + tt~ angular distribution. The -k + tt^ amplitude 
contains a broad pole at (440 — 470) MeV. In some sets of data it appears as a peak, e.g. in 
J/\& — > ui7T + 7r~ [18J. In others it appears as in elastic scattering. In principle, it therefore needs 
to be fitted by an amplitude [Ai + A 2 (s — sa)]/D(s), where sa is the Adler zero, 0.0097 GeV 2 



4 



and A are coupling constants. An algebraic parametrisation of the denominator D(s) is given 
in Ref. [19J. It is not correct to fit the data with the E791 parametrisation, a Breit-Wigner 
amplitude of constant width. That has a phase variation in complete disagreement with tttt 
elastic scattering and violates Watson's theorem |20j. Furthermore, it contains an unphysical 
pole just below the irn threshold. 

3 X,Y,Z 

It is well known to theorists that the opening of strong thresholds can move resonances by large 
amounts. As one example, Barnes [21] has studied charmonium loop diagrams and concluded 
that mass shifts are potentially > 100 MeV; he finds that they tend to be similar for all states 
in a multiplet with similar radial excitation number n but different orbital angular momentum 
L, e.g. the n = 2 3 -Po,i,2 an d 1 P\ states. 

The Z(3930) is identified clearly in Belle data for e + e~ — > e + e~ + DD by a strong P-wave 
component in DD with mass M = 3929 ±5 ±2 MeV, T = 29±10±2 MeV. [22]. It is now adopted 
by the PDG as the x C 2(2P) radial excitation of Xc2 P3]. There are two observations still labelled 
X(3945) by the PDG. One is by Babar in B -> K(uJ/^) with M = 3914.6±|;f ± 2.0 MeV, 
r = 34lg 2 ±5 MeV [23]. The second, with less events (49±15), is by Belle in e+e" e + e-(wJ/f) 
with M = 3916 ± 3 ± 2 MeV, T = 17 ± 10 ± 3 MeV [21]. The mass difference from Z(3930) is 
14 ±6 ±4 MeV, so it is quite possible this is the Xc2(2P) decaying to uJ/ty. Another possibility 
is the Xco(2P). 

There is considerable scatter in masses predicted for these states. A representative collection 
of recent predictions is shown in Table 1. Mass differences between 3 P 2 and 3 Po radial excitations 
are large, but can be affected by thresholds and form factors. If the n = 2 3 P 2 predictions can 
be trusted, at least approximately, this state has been pulled down in mass by a large amount. 
The strong signal for uJ at 3915 MeV is a hint of mixing between cc and the ojJ/'fy channel 
and attraction to this threshold and/or the DD* threshold. 



State 


BGS [25] 


EFG [26] 


[27] 


3 P 2 


3972 


3972 


3941 


3 P 


3925 


3929 


3900 


3 P> 


3852 


3854 


3839 


l Pi 


3934 


3945 


3909 


X D 2 


3799 


3811 


3799 



Table 1: Recent predictions of masses of n = 2 3 P states and n = 1 1 D 2 in MeV. 

The V(3945) claimed by Belle [28] at 3943 ± 11 ± 13 MeV with T = 87 ± 22 ± 26 MeV does 
not agree well with masses near 3915 MeV. However, they are in the same decay channel and 
with a small stretch of the errors may both be consistent with Z(3930). 

A puzzle is where the 3 Po state is. Its obvious decay mode is to DD. Its width is strongly 
dependent on the form factor for the DD channel [29], [30]. An early paper of Babar comments 
that separation of DD events is obscured by mis- identified DD* events, so it could have escaped 
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detection just below that threshold [31]. The other obvious possibility is that it could be the 
3915 MeV peak in ujj/ty, but then the inference is that it couples strongly to that threshold. 

The X(3940) is observed by Belle in e + e~ — > J/^J + an inclusive collection of DD and DD* 
events. From 52 identified DD* events, Belle find a mass of M = 3943i£ ±6 MeV [32]. There arc 
two obvious possibilities for its spin-parity. Firstly, it could be the missing 1 Pi state. Secondly, 
it could be Xc2(2-P) decaying to DD* with orbital angular momentum I = 2, though there is 
presently no claim to observe it in DD. The rising phase space for i = 2 would naturally 
move the peak up in mass. To illustrate this, suppose Z(3930) is described by a Breit-Wigner 
resonance of constant width with Belle parameters. A simple calculation using DD* phase space 
and an i = 2 centrifugal barrier for decay with a radius of interaction of 0.73 fm [33J gives a 
peak at 3939 MeV with width 43 MeV; this is close to Belle parameters, but perhaps fortuitous 
in view of experimental errors for both Z(3930) and X(3940). Also strong coupling to the DD* 
channel would move the peak by an amount which cannot be calculated without knowledge 
of the relative branching ratios to DD and DD*. A puzzle is that Belle do not observe any 
significant B ->■ K(DD*) signal near 3940 MeV [33]. 

The spin-parity analysis of e + e _ — > J/^ + 2 + involves spins 1~ — > l _ + 2 + + L; since momenta 
in the final state are large, L can run up to 4. This makes a complete partial wave analysis 
impossible with present data and J p = + can be confused in this reaction with the helicity 
component of the 3 P2 state. 

4 Spin-Parity Analysis 

Spin-parity determinations of these states are needed if any progress is to be made in under- 
standing them. The cleanest route to J p determination lies in B — > K + X, X — > uJ/^fr (or 
pJ/'ty) or DD*. This is because both B and K are spinless, with the result that the kaon 
is produced with known orbital angular momentum equal to the spin of X; this simplifies the 
partial wave analysis greatly. However, that analysis must use formulae which describe both 
production and decay of X and must use spin information from both uj and J/vP or from the 
D* . The spin of the u is normal to its decay plane in the rest frame of 7r + 7T~7r°. The J/ip is 
detected via decays to lepton-lepton, and its polarisation vector e is normal to the lepton-lepton 
axis in its rest frame. These are key pieces of information. With them, matrix elements for the 
full amplitude contain correlations between the spectator kaon and decays of u and J/^> which 
are highly distinctive; without these correlations, most of the relevant information is discarded. 
Formulae for these matrix elements were published in 2005 [H] and will be reviewed here in 
outline, together with comments on the simplest choice of reference frames. 

If formulae are written in terms of vectors and tensors, the vital features are transparent. 
The formulae can be written in any reference frame, but in practice they simplify greatly in the 
rest frame of X,Y,Z for reasons explained algebraically in Appendix 1. A general point is that 
all matrix elements are orthogonal; that is a valuable check on computer programmes. 

The decay of the resonance X is described by combined spins s of u and J/ty (or D + D*). 
The Z(3930) is likely to decay dominantly with i = 0, though I = 2 is allowed in principle. 
However, the centrifugal barrier factor for t = 2 is quite strong. Table 2 exhibits Blatt-Weisskopf 
centrifugal barrier factors for the amplitude using a generous radius of interaction 0.73 fm, 
including the convolution of both particles in the final state; (average masses are used for D 
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IVlcLoo I IVlt: V J 


7~) 7~)* rlQT'T'lfiT' 

LJ LJ Udlllt;! 




3900 


0.1905 


0.1105 


3915 


0.2880 


0.1835 


3930 


0.3679 


0.2498 


3945 


0.4361 


0.3096 


3960 


0.4937 


0.3635 



Table 2: L = 2 centrifugal barriers for DD* and uJ/ty. 



and .D*). For a smaller radius, the barrier factors are even smaller. The I = 4 barrier makes 
that amplitude negligible. The primary objective of the analysis should be to determine the 
dominant amplitudes, hence should start with J p = + and 2 + , both with £ = 0. Further £ = 2 
amplitudes may then be introduced one by one with penalty functions which constrain their 
magnitudes unless the data really demand them. This may be done, for example, by adjusting 
the penalty function to contribute \ 2 or l°g likelihood in the range 1-4. 

Y 

x< T 




Figure 3: Axes describing B — > KuJ/ty. 



A simple choice of axes is shown in Fig. 3 with the Z-axis along the direction of the J/\l/ and 
Y along the normal to the plane of K, u and J/\l/. The formula for u> — > 7r + 7r~7r is a 4-vector 
W given in Lorentz invariant form in Section 4.2 of [H]. This formula is easily evaluated in any 
reference frame, such as that of Fig. 3. It includes naturally the well-known enhancement near 
the edge of the Dalitz plot for its decay to 7r + 7r~7r°. 

The J/ty decay is simplest in its rest frame. The effect of the Lorentz transformation to the 
rest frame of X,Y,Z is small because the J/\l/ is heavy and therefore non-relativistic in that 
frame. The gymnastics required for the Lorentz transformation of its polarisation vector e are 
given in Appendix 2. 

For J p = 1 + , the kaon in the production process carries L = 1. In the resonance rest 
frame, this angular momentum is described by the 3-momentum K of the kaon. In the decay 
X(3872) — > pJ/ty, the spin 1 of the p is represented by the vector P = k(TTi) — A;(7r 2 ), where k 
are momenta in any reference frame. In the non-relativistic limit the full matrix element reduces 
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to M = e.P A K . There is a simple trick which can be used to do the spin average over e. This 
is to write e x = 1 and e y = i = \/— 1; then intensities are obtained by taking the modulus 
squared of M. For X(3872), £ = 2 in the decay is eliminated by the tiny width of the state. 
The threefold dependence on the lepton axis, the spin of the p and the kaon direction is then 
unique and highly distinctive. 

For J p = + , the matrix element is 

M = W.e = W 1 e 1 + W 2 e 2 + W 3 e 2 - W 4 e 4 , (6) 

where W and e are 4-vectors. Without the last term, which is very small, the matrix element 
is just the scalar product of 3- vectors W and e. With this simplification, the intensity is sin 2 a, 
where a is the angle between W and the lepton axis. The simplicity of this formula may actually 
make the + state harder to identify than other amplitudes. 

For J p = 2 + , the kaon carries L = 2 and is described by a tensor 

T a p = K a K f) -{l/3)6 a f ) (K*). (7) 

Here 5 is the usual Kronecker delta function and K 2 = K\ + K\ + K 2 — K\. For decays with 
I = 0, the combined spin of u and J/ty is s — 2. The decay is then described by the tensor 

T a p = e a Wp + £f) W a - (2/3) (e^) (8) 

in the resonance rest frame; e M W implies the summation over /i = 1 to 4 and takes the same 
general form as K 2 above. The full matrix element is just the contraction T a pT a P. It is very 
distinctive because of the dependence on angles of all of K, W and e. 

A simulation of these J p , assuming full geometrical acceptance, suggests that 50 events would 
give a separation at 90% confidence level between + , 1 + and 2 + if only one of them contributes 
a significant amplitude. Babar and Belle between them have many more events than this, with 
the exception of X(3940) — > DD* . If more than one J p contributes, the separation is obviously 
poorer, and an overlap between + and 2 + could be a difficult case. Formulae for i = 2 in decays 
are given in Ref. [14] . 

A distinctive feature of X(3872) is that its mass lies within 0.3 MeV of the D q Dq threshold. 
There is a simple explanation. It is attracted to this threshold by the dispersive effect due to 
the opening of the threshold. Lee et al. [35J model the A(3872) including exchange of 7r and 
a between D and D*. They solve the Schrodinger equation, including the DD* channel; this 
is equivalent to evaluating the dispersion integral of Eq. (2), since their solution is explicitly 
analytic. One should not be misled into thinking that the effect of DD* is small because few 
events are seen in X(3872) — > D Dq] there are few events because of the narrow width of 
X(3872). However, the attraction to the threshold does not involve the resonance denominator 
in Eq. (2); it arises from the integral over all phase space for DD*. A parallel example is the 
deuteron. Its binding energy is not governed by its tiny width; it comes from virtual meson 
exchanges over the whole left-hand cut. 

Present data have failed to distinguish between J p = 1 + and 2~ for A (3872). The key 
point is to check whether partial wave analysis is consistent with J p = 1 + or not. If it is, the 
fact that the mass coincides with the DD* threshold is significant information. A state with 
J p = 2~ couples to DD* with I = 1 and there is no obvious reason why it should lie close to this 
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threshold. In fact, predictions for the mass of 77^(1 1 D2) in Table 1 lie far below X(3872); this 
is because spin-splitting for D states is smaller than for P-states and the mass of the ?/>(3770) 
(mostly n = 1 3 Pi) is well known. For J p = 2~, formulae are more messy and are given in Ref. 
P3]. There are are in principle three amplitudes, but one requires I = 3, which is irrelevant for 
the narrow X(3872). The others have I = 1 and combined spins either s = 1 or 2. 

A full partial wave analysis of Belle data for X(3940) — > DD* will in general be ambiguous 
with only ~ 50 events. However, this may be enough events to test the possibility that J p = 0~. 
The spin information is carried by the vector 

Q = k D -K- ° n [k D + K], (9) 

M D * 

where k are momenta (in any frame). With J p = 0~, production requires L = 1 and decays 
to DD* with I = 1, so only a single amplitude is present; it has the distinctive feature of 
going through zero at the centre of the angular distribution for decays. There is a missing 
1 Pi cc state. From predicted masses in Table 1 and the fact that the mass of the n = 1 1 Pi 
state is slightly about the 3 Pi state, it seems likely that the x Pi state will lie above X(3872). 
However, it is not easy to find. It has charge conjugation C = —1, while ojj/'fy has C = +1, 
so this decay is forbidden. It can appear in DD*. However, it appears in that channel in the 
combination D + D*_ + D Dq (because of the minus sign associated with D ), while 3 Pi decays 
with the opposite relative sign. Without interferences with a state of established J PG , 1 P± and 
3 Pi are not separated in any particular charge state. 

5 Conclusions 

It is necessary to be specific about what tetraquark configuration is meant in every discussion 
of them. Meson-meson combinations mixed with qq are to be expected via t- and w-channel 
exchanges and also from the s-channel decay. 

The increase in branching fractions for e + e~ — > tt + tt~T(1S,2S,3S) as the mass increases 
from T(45*) to T(10860) may be understood qualitatively via the opening of real decays to 
BB*, B*B* and B S B* channels in the mass range 10559 to 10826 MeV, followed by de- 
excitation to 7T7rT. Observed branching fractions are not seriously out of line with expec- 
tation. The composition (bn)(bn) may be tested by measuring the intensity ratio given by 
a[T(lS)K + K-]/a[T(lS)K°K°), which Ali, Hambrock and Mishima show to be 1/4 for exotic 
tetraquarks; for a BB* composition or bb, this ratio should be 1. 

Spin-parity analyses of X,Y,Z states are needed and can be done straightforwardly in several 
cases using formulae developed in Ref. [H]. It is suggested here that they can be explained 
naturally as the expected n = 2 radial excitations of 3 Pi, 3 P2 and perhaps 3 Pq. It is essential to 
make use of the full forms of matrix elements for production and decay for partial wave analysis. 
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7 Appendix 1 



Chung gives an introduction to tensor formalism [36]. Foundation work was done by Zemach 
[37] . The convention used here is that the Lorentz metric g^ v is such that four- vectors are 
written as (—pi, —p 2 , ~P3, E = p±), i.e. g^ v has only diagonal elements (-1,-1,-1,1). 

The decay of a scalar particle, e.g. o — > mr, is described by its mass {k\ + ^2)4 = -P4; here k\ j2 
are 4-vectors of the decay pions. A spin 1 particle needs to be described by a vector constructed 
to be orthogonal to the scalar expression. This is done, for example in describing p — » ittt, using 

(fcuV = aivih + k 2 ) v , (10) 

where 

9^ = 9^-{PM/Mi- (11) 

Substituting (11) into (10), 

Vti, = a-ft - h). - Pl " hUh + J? )A(kl + k2) " - (12) 

p 

The orthogonality between vector and scalar is demonstrated from the dot product of (k^)^ 
with (ki + k 2 ) fJ/ : 

g^(ki - hUh + h)^ - {Ml - M 2 2 )M 2 p /M 2 p = 0. (13) 

In the rest frame of the decaying particle, g^ u has diagonal elements [-1,-1,-1, l]-[0, 0,0,1] = 
-[1,1,1,0]. This gives the important result that k is a 3- vector in the rest frame of a resonance. 
Then k describes the angular momentum carried by the particle. Tensor expressions for higher 
spins are built up as described by Zemach and Chung so that they are orthogonal to lower spins. 
They are symmetric and traceless, i.e. diagonal elements add to 0. We shall need only the tensor 
T^y for spin 2. 
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Suppose the lepton pair from J/^ decay has components (sin 9 cos 0, sin 6 sin 0, cost 
in the J/\l/ rest frame parallel to X, Y, Z of Fig. 3. This vector may be expressed as 



/ COS ( 

sin ( 

V o 




cos 9 sin 9 

1 
— sin 9 cos 9 




m axes 



(14) 



The spin e of the J/^f is orthogonal to this lepton axis. In axes x,y,z aligned along the lepton 
axis, e is given by a vector (cos R, sin R, 0) and eventually it is necessary to average over R. In 
X,Y,X, axes, e is given by the converse relation to Eq. (14) with the result 



/ cos 9 
1 
^ sm9 




/ COS 

— sin ( 

V o 




(15) 
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cos 6* (cos cos R + sin sin R) 

— sin cos i? + cos sin i? |. (16) 
sin #(cos cos R + sin sin R) 

Under the Lorentz transformation back to the centre of mass frame, this becomes 



/ cos 6 (cos cos R + sin sin R) \ 
— sin cos i? + cos sin i? 
7 sin #(cos cos i? + sin sin i?) 
V /37sin#(cos0cos.R + sin sin i?) / 



(17) 



where 7 and /3 are the usual expressions for the Lorentz transformation. This approach is simpler 
than that given in Ref. [H]. 
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